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We consider three fluid phases meeting at a line of common contact and study the linear excesses per unit 
length of the contact line (the linear adsorptions Aj) of the fluid's components. In any plane perpendicular 
to the contact line, the locus of choices for the otherwise arbitrary location of that line that makes one of the 
linear adsorptions, say A2, vanish, is a rectangular hyperbola. Two of the adsorptions, A2 and A3, then both 
vanish when the contact line is chosen to pass through any of the intersections of the two corresponding 
hyperbolas A2 = and A3 = 0. There may be two or four such real intersections. It is required, and is 
confirmed by numerical examples, that a certain expression containing the adsorption of component 

1 in a frame of reference in which the adsorptions A2 and A3 are both 0, is independent of which of the 
two or four intersections of A2 = and A3 = is chosen for the location of the contact line. That is 
not true of A 1 ( 2 ,3) by itself; while the adsorptions and the line tension together satisfy a linear analog of 
the Gibbs adsorption equation, there are additional, not previously anticipated terms in the relation that are 
required by the line tension's invariance to the arbitrary choice of location of the contact line. The presence 
of the additional terms is confirmed and their origin clarified in a mean-field density-functional model. The 
additional terms vanish at a wetting transition, where one of the contact angles goes to 0. 
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1 Introduction 



It was long ago remarked by Buff [1] (as quoted by Melrose [2] and Good [3]), and Lovett [4], that a mea- 
surement of the temperature or chemical-potential dependence of the interfacial tension between two phases 
can yield a length characteristic of the microscopic structure of the interface; viz., the distance between two 
well-defined Gibbs dividing surfaces. It is then a microscopic length obtainable from a macroscopic mea- 
surement, which is remarkable. What makes it possible is that the measurement is made in three dimensions 
while what is measured is a property of the interface, and thus of a system of one less dimension. The result 
is then the characteristic length of the "missing" dimension, which is microscopic. It is a similar reduction 
in dimension from three to two that made possible Rayleigh's famous estimate of the length of a molecule 
by measurement of the area of a monolayer that had spread from a known volume [5]. 

The purpose of the present paper is to give a geometric interpretation of the temperature or chemical- 
potential dependence of the tension of the line of contact of three phases, analogous to that found by Buff 
and Lovett for the interface between two phases. Here we first rederive the result of Buff and Lovett, and 
we then state the analogous question for a three-phase contact line. 

Let the tension of the interface between two phases a and (3 be <T Q/ g, the chemical potentials of the various 
chemical components i of the system be pi, and the densities of those components in the bulk phases be pf 
and pf . It is convenient [6], for simplicity of notation, to include the temperature T among the pi and the 
entropy densities s a and s 13 among the pf and pf , so that in a c-component system the index i = 1,2,..., 
c+1. 

In Fig. 1 is shown, schematically, the two phases a and (3 with a diffuse interface between them, and an 
arbitrary Gibbs dividing surface chosen to be at the height zq on the vertical scale of distance, z. Let Tj be 
the adsorption (surface excess per unit area) of the component i with respect to that dividing surface. Then 
the Gibbs adsorption equation is 

i 

Of the c + 1 chemical potentials pi (among which is the temperature T) only c are independent because they 
are related by the Clapeyron equation [6] 

£ (P?-P?)dp t = 0. (1.2) 

i 

The adsorption Ti(z) depends on the location zq of the dividing surface via 

A T, = (pf - pf ) Az (1.3) 

where AT{ is the change in T{ when the dividing surface is displaced by the distance Az in the direction 
from phase [3 towards phase a [1-4,6]. Let Zi and Zj be the locations for the dividing surface that make 
Ti = and Tj = 0, respectively, and, in standard notation, let be the adsorption Tj when the dividing 
surface is chosen to be at Zj, where Tj = 0. Then from (1.3), 

%)= (pf-p f l)(z 3 -z l ). (1.4) 
In a one-component system, or in a c-component system at fixed p3, . . . , p c +i, there is only one inde- 
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pendently variable ^, say fi±. Then from (1.1) and (1.4), 



-T^- = (Pi - Pi) (zi - ^2) (fixed /x 3 , . . . , // c +i); (1-5) 
"Ml 

i.e., from a measurement of da a p/dfii and the known densities in the bulk phases one obtains both the 
magnitude and sign of the distance z\ — z 2 between the two dividing surfaces of T\ = and r 2 = 0. This 
is the Buff-Lovett principle [1-4,6]. 

In Fig. 2 we show three phases a, (3, 7 meeting in a line of common contact, which is where the three 
interfaces a/3, f3j, and 07 meet. We do not consider curvature effects, so the interfaces are planar and the 
contact line is straight. The dihedral angles between interfaces are the contact angles, which we also call 
a, [3, and 7, with a + (3 + 7 = 2ir. The location, although not the direction, of the contact line is arbitrary, 
just as is the location of the dividing surface in Fig. 1. In Fig. 3 is a view along the contact line, which is 
perpendicular to the plane of the figure, so it appears as a point, the interfaces appearing as lines meeting in 
that point. Two alternative and equally arbitrary choices for the location of the contact line, and so also for 
the locations of the interfaces, are shown in the figure. There is thus a two-dimensional infinity of choices for 
the location of the contact line; i.e., it may be chosen to pass through any point in the plane of the figure. In 
the analogous Fig. 1 for the interface between two phases there is only a one-dimensional infinity of choices 
of dividing surface. The line tension r, like the surface tensions o a p, etc., is independent of what choice is 
made. 

There are linear adsorptions Aj (linear excesses of the components % per unit length of the contact line) 
analogous to the surface adsorptions Tj, and there are curves in the plane of Fig. 3 that are the loci of the 
choices for the location of the contact line for which any of the linear adsorptions Aj is 0. We ask what those 
loci are, and what geometrical information is contained in A 1 ( 2 ,3)» the linear adsorption of component 1 in 
a frame of reference in which A2 = A3 = 0; i.e., when the contact line is chosen to be at the intersection 
of the curves A2 = and A3 = 0. Likewise, we ask for the geometric interpretation of dr/dfii, and of the 
connection between dr/dfi\ and A 1 ( 2 ,3). The answers will be the analogs for the three-phase contact line of 
the relations (1.1), (1.4), and (1.5) for two-phase interfaces, but will prove to be a more complicated story. 

In Section 2 we find that the locus of choices for the contact line that make Aj = (briefly, the curve 
Aj = 0) is a rectangular hyperbola, and that any two such hyberbolas, say A2 = and A3 = 0, may intersect 
in two or four points. We find, however, that A 1 ( 2 ,3) depends on which of the intersections is chosen as the 
location of the contact line. This means that — A 1 ( 2 ,3) cannot be drjd^\, which is invariant to the arbitrary 
location of the contact line; i.e., that the linear analog of the Gibbs adsorption equation (1.1) must contain 
hitherto unrecognized terms that restore the invariance. This is the subject of Section 3. The additional 
terms are found to vanish at a wetting transition, which is where one of the contact angles closes down to 0. 
In Section 4 we confirm in a simple mean-field density-functional description what in Section 3 we suggest 
is the physical origin of the extra terms, and we see how in principle these terms would be calculable from 
a microscopic theory. The results are briefly summarized in the concluding Section 5. 



2 Linear adsorption Aj, locus Aj = 0, and K(j,k) 

In Fig. 4 we show the Neumann triangle associated with the three-phase equilibrium. This is in the plane of 
Fig. 3, with a particular choice of location of the contact line. The sides of the triangle are perpendicular to 
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the interfaces a(3, (3j, and cry, intersecting them at the respective distances R a p, Rpy, and R a y from the 
contact line. The vertex angles of the triangle are the supplements it — a etc. of the contact angles. The 
three interfaces divide the triangle into the three areas A a , Ap, Ay, as shown. 

We introduce the unit vectors e a p, ep 7 , e a7 in the directions of the interfaces, and another set of unit 
vectors g Q/ 3 , g/3 7 , g Q7 that are respectively perpendicular to the e a p, etc., and are rotated from them coun- 
terclockwise. These are illustrated in Fig. 5. There will later also occur dyads of the form e a p g a p, etc. 

With the Neumann triangle fixed in space, the distances i? Q| g(r) etc., areas A a (r) etc., and surface 
adsorptions T" /3 (r) etc., are all functions of the arbitrary location r of the contact line. For two different 
choices ro and r they are such that 

R a p(r) - R al3 {ro) = -e a p ■ (r - r ), (2.1) 

and similarly for Rp 1 and i? Q7 ; 

A a {r) - A a (r ) = - R a/3 {ro) g a p ■ (r - r ) + i? a7 (r ) g a y ■ (r - r ) 

+ \( r ~ r o) • (e Q( 9 ga/3 - e a7 g Q7 ) • (r - r ), (2.2) 
and similarly for Ap and Ay with a, j3, 7 permuted cyclically; and, from (1.3), 

If (r) - Ff (r ) = (pf - pf ) g a p ■ (r - r ), (2.3) 

and similarly for rf 7 and T" 7 with a, j3, 7 permuted cyclically. Note that these are invariant to the inter- 
change of r and ro- That is obvious for (2.1) and (2.3), and for (2.2) it is a consequence of (2.1). Note 
also from (2.2) and its analogs for Ap and Ay obtained from A a by cyclic permutation of a, j3, 7 that 
A a + Ap + Ay is independent of the location of the contact line, as it must be because it is just the total area 
of the Neumann triangle. 



It is to be appreciated that while the contact line may be located arbitrarily at any r in the plane of Fig. 4, 
the properties of the system are spatially varying and, with the physical system itself fixed in space, different 
choices of r correspond to physically different structures about the contact line. While T" /3 (r) does vary 
with r, it is a well defined and in principle measurable function of it. 



That is true also of the linear adsorptions Aj, which are functions A, (r) of the choice of location r of the 
contact line with respect to which they are defined. For any two such choices ro and r, 

A,(r) - A,(r ) = - {pf [A a (r) - A a (r )} + • • • 

+ [Ff (r) R ap (r) - Tf (r ) R a p(r )] +■■■} 

where the omitted terms marked by • • • in the first line have (3 and 7 in place of a and those in the second 
line have (3j and 07 in place of a(3. From (2.1)-(2.4), then, 

Ai(r) - A,(r ) = [rf (r )e Q/3 + rf 7 (r ) e/37 + T^(r )e a y] ■ (r - r ) 

+ - r o) • [(p? ~ Pi)e a p g a p + (Pi - p])epy gpy (2.5) 

+ (Pi ~ Pi)eay ga 7 ] " ( r ~ r o)- 
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Note that there is no remaining reference to the distances R a p etc.: the difference Aj(r) — Aj(r ) depends 
only on the locations ro and r in the physical system and makes no reference to the Neumann triangle. Note 
also that (2.5) is invariant to the interchange or r and ro, as a consequence of (2.3). 

When the "component" i is the grand-canonical free energy, its linear adsorption A, is the line tension r, 
its surface adsorptions etc. are the interfacial tensions a a/ 3 etc., and its bulk-phase densities pf , pf , p] 
are all equal to —p, the negative of the uniform bulk-phase pressure (since the interfaces are assumed planar 
and the contact line straight). Then (2.5) together with the condition 

e Qj a + crpj ep 7 + o ai e Q7 = (2.6) 

of mechanical equilibrium yields r(r) — t(tq) = 0, which is the required invariance of the line tension to 
the location of the contact line. 

If we have a system of rectangular coordinates x, y with origin at ro and with the x and y axes in the 
directions of e a7 and g a7 respectively, the equation Aj(r) = 0, from (2.5), is the conic 

A iV 2 + B iX y + dx 2 + D iV + E iX + Fi = (2.7) 

with 

-d = A = -l(pf - p<l ) sin 2a + - pj) sin2 7 

Bi = \{pl~ Pf) + \{P?~ P?) «* 2a + i(pf - pj) cos 2 7 

A = -rf (r ) sin « + rf (r ) sin 7 

Ei = Tf (r ) cos a + rf(r ) cos 7 + r^(r ) 

F^AiCro), (2.8) 

where a and 7 are two of the contact angles: e a p ■ e ai = cos a and e^ 7 • e a7 = cos 7. Because Cj = —Ai 
the conic is a rectangular hyperbola (perpendicular asymptotes). Although the form of its equation depends 
on the choice of origin ro and the chosen orientation of the x and y axes, the hyperbola itself does not; it 
is a physically well-defined curve with a well-defined location and orientation in any plane perpendicular to 
the direction of the contact line, and could in principle (although perhaps only with difficulty in practice) be 
traced by measurements on the fixed equilibrium system. 

Let Aj(r) = be similarly defined as the locus of those choices r for the location of the contact line 
that make the linear adsorption Aj = 0. The two hyperbolas Aj(r) = and Aj(r) = may in general 
make two or four real intersections with each other. They cannot make more than four, as can be seen by 
translating and rotating axes so that one of the hyperbolas becomes of the form xy = constant in the new 
coordinate system. Then substituting y = constant jx in the equation of the other and multiplying through 
by x 2 yields a polynomial equation in x of the fourth degree for the locations of their intersections, which 
can have at most four real roots. No intersections, or confluent intersections, are possible but non-generic. 

On dividing (2.7) through by Ei and recalling that Cj = —Ai, we see that the hyperbola Aj(r) = is 
determined by four parameters, A{/Fi, Bi/Fi, Di/Fi, and Ei/Fi. Thus, typically only one such hyperbola, 
not two distinct ones, could pass through four given points, but two distinct ones may do so if the four points 
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in question, (x\, y±), . . . , (x^, 2/4), are related by 
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for then the four points would no longer determine the four parameters A,jFi etc. uniquely. Thus, either 
the two rectangular hyperbolas Aj(r) = and Aj(r) = intersect each other in two points, or in four, and 
in the latter case the four intersections satisfy (2.9). 

Figure 6 illustrates the case of two intersections of the hyperbolas A2 = and A3 = 0, with the 
parameters 

A 2 /F 2 = -1, B2/F2 = 9/10, D2/F2 = E2/F2 = 3 

^3/^3 = 0, B 3 /F 3 = -l/5, D 3 /F 3 = E 3 /F 3 = 0. (2.10) 

Figure 7 illustrates the case of four intersections, with 

A2/F2 = -1, B2/F2 = D2/F2 = E2/F2 = 
A3/F3 = -1/10, B 3 /F 3 = -l/5, D 3 /F 3 = 2/5, E 3 /F 3 = -l. (2.11) 

Values of the parameters on the right-hand sides of (2.8) that are implied by these values of Ai/Fi etc. are in 
the Appendix, as are numerically accurate values for the x, y coordinates of the various intersections. The 
coordinates of the four intersections in Fig. 7 satisfy (2.9). 

Now let Yjk be any one of the two or four intersections of Aj(r) = with Afc(r) = and let r-j be any 
point on the locus A«(r) = 0. Then from (2.5), taking r to be r-j and 1*0 to be rjk, we obtain for Aj^), the 
linear adsorption Aj in a frame of reference in which the contact line is chosen to be at a point where the 
linear adsorptions Aj and A& are both 0, 

A i(j,k) = ~ [rf (r jk )e aP + rf 7 (r jfc ) e/37 + T^(r jk )e aj ] ■ (r t - r jk ) 

~ ~ r jk) ■ [(p? ~ Pi)e a p g«/3 + (pf - pJ)e Pl g/3 7 ( 2 - 12) 

+ (p] - Pi)e aj ga-y] ■ (*i ~ *jk)- 

This can be written in several alternative ways by taking advantage of the invariance of (2.5) to the inter- 
change of r and ro, but the form in (2.12) is as convenient as any. 

Equation (2.5), or the nearly equivalent (2.12), is fundamental in our analysis. Equation (2.5) is the 
analog for the contact line of (1.3) for an interface and (2.12) is the analog of (1.4). In the following section 
we consider the linear analogs for the line tension of the Gibbs adsorption equation (1.1) and the Buff-Lovett 
principle (1.5), and the roles of Aj and A i(j fc ) in these relations. 
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3 Gibbs adsorption equation for the contact line 



One reads in Gibbs [7], "These [three-phase contact] lines might be treated in a manner entirely analogous 
to that in which we have treated surfaces of discontinuity. We might recognize linear densities of energy, of 
entropy, and of the several substances which occur about the line, also a certain linear tension. With respect 
to these quantities and the temperature and potentials, relations would hold analogous to those which have 
been demonstrated for surfaces of discontinuity." By following this prescription, exactly analogously to the 
derivation of the adsorption equation (1.1), the linear adsorption equation in the form 

dr = - ^ At dm (3.1) 

i 

was found in earlier work [8,9]. We shall now see that, unlike the surface adsorption equation (1.1) for da, 
this is not invariant to the choice of location of the contact line, as it must be, and as was verified in the 
remark below Eq.(2.6). Something is missing from (3.1). 

In three-phase equilibrium there is in addition to (1.2) a second Clapeyron equation [8], 

£ {fJ 3 i -pl)dm = Q, (3.2) 

i 

where, as in Section 1, the entropy densities s a etc. in the bulk phases and the temperature T are included 
among the pf etc. and the pi, respectively. Because of (1.2) and (3.2) together, there are now only c — 1 
independent p, L . 

With the Gibbs adsorption equation (1.1), which holds whatever the dividing surface with respect to 
which the T L are defined, and the Clapeyron equations (1.2) and (3.2), it follows from (2.5) that with in- 
finitesimal changes dpi, 



^2 Aj(r) dpi - ^2 Aj(r ) dpi = - (e Q/3 da af3 + da Pl + e ay da ay ) • (r - r ). (3.3) 

i i 

Only if e Q/ 3 da a p + e/3^ dapy + e a7 da a7 were identically would the left-hand side of (3.3) be for all 
r and ro, so only then would £ Aj dpi be independent of the location of the contact line, as dr is required 
to be. Equivalently, from (2.6), this would require a a p de a p + de^ + a ai de ai to be identically 0. 
We may arbitrarily hold the direction of one of the interfaces, say e Q7 , fixed while making the infinitesimal 
changes dpi. Then in terms of the two orthogonal unit vectors e a7 and g a7 and the contact angles a and 7, 

<r a /3 de a/3 + a de^ + a ai de ai = (a^ cos 7 dj - a a/3 cos a da) g a7 ^ 

— (crg 7 sin 7 c?7 + a a p sin a da) e Q7 . 

For this to vanish would require the coefficients of both g Q7 and e Q7 to vanish; and so, with a a f3, a^, and 
da/dj all generally finite and non-zero, it would require that cos a sin 7+sin a cos 7 = 0; or sin(a+7) = 0; 
or sin/3 = with (3 = 2tt — (a + 7) the third contact angle. Thus, the left-hand side of (3.3) will vanish 
identically only when (3 = 0; i.e., only when the (3 phase wets the interface between the a and 7 phases (or 
when any of the phases wets the interface between the other two, (3 having been singled out only by having 
imagined e Q7 fixed). 
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The non-invariance of E Aj dfii may be illustrated with the intersecting hyperbolas in Figs. 6 and 7. Let 
H4, . . . , /i c +i be fixed, so that, with the two Clapeyron equations, there is only one independently variable 
Hi, say fi\ again, as in Section 1. Then E Aj d/Xj = A 1 ( 2j 3) d^i when A 2 = A3 = 0, so the non-invariance 
of E Aj d/Xj would be illustrated by the non-invariance of ^1(2,3) as calculated from (2.12); i.e., by the 
dependence of A 1 ( 2j 3) on just which of the two or four intersections of A 2 = with A3 = in Fig. 6 or 7 is 
identified as r 2 3. For the purposes of this illustration we take for the hyperbola Ai = the parameters 

A 1 /F 1 = -2/5, B 1 /F 1 = -5/2, D 1 /F 1 = -7, E 1 /F 1 = -l. (3.5) 

In Fig. 8 we show this hyperbola along with the two points PI and P2 from Fig. 6, and in Fig. 9 we show 
the same hyperbola (3.5) along with the four points PI, . . . ,P4 from Fig. 7. 

Values of the physical parameters p",r°^(ro), etc. in (2.8) that go with the Ai/Fi etc. in (2.10), 
(2.11), and (3.5) and are thus required for the evaluation of ^1(2,3) from (2.12), are in the Appendix, along 
with numerically accurate values of the coordinates x, y of the intersections PI and P2 in Figs. 6 and 8 
and of PI, . . . , PA in Figs. 7 and 9. When r 2 3 in (2.12) is taken as the point PI in Figs. 6 and 8 one 
finds A X ( 2 3 ) = 5.90135 for every point ri on both branches of the hyperbola Ai = 0, whereas when r 2 3 is 
identified with P2 one finds the different value ^1(2,3) = —56.4125 for every point ri on Ai = 0. Similarly, 
when r 2 3 is identified in turn with the four points PI, P2, P3, P4 in Figs. 7 and 9 one finds for ^-1(2,3) the 
respective values -7.17944, -15.694, 9.73286, 136.541, each of these for every 1*1 on both branches of 
Ai = 0, but differing from each other. The non-invariance of ^-1(2,3) in these examples, illustrating that of 
E Aj d\ii, is thus manifest. 

What is missing from (3.1) are terms that effectively cancel the non-invariant parts of E Aj d[ii. In their 
general treatment of the thermodynamics of line tension, Boruvka and Neumann [10] introduced reversible 
work terms involving changes of the contact angles (as well as curvature terms, which we have been sys- 
tematically omitting). We now include such terms in the infinitesimal energy change dU, but since they are 
absent from the bulk and surface energies they emerge only in dll e , the infinitesimal change in the linear 
excess energy. For length L of the three-phase contact line, we thus take dU e to be 

dll e = ^2 A*i dN i +r dL + L(c a da + cp d(5 + c 7 dry), (3.6) 

i 

where N% is the linear excess of component i and where, following our convention, the linear excess entropy 

S £ is included among the Nf with T the associated /ij. It is the terms L(c a da-\ ) that were absent in the 

earlier treatments [8,9]. Since a + (3 + 7 = 2ir, only two independent combinations of the three coefficients 
c a ,cp, c 7 occur; e.g., c a — cp and eg — c 7 . 

These coefficients depend on the location of the contact line. As remarked earlier, different choices put 
the contact line in physically different environments in the equilibrium system. Therefore the same contact 
angle changes da, d(3, dj would produce different distortions in the spatial variations of the component 
densities, and so induce different changes in dU , for different locations of the contact line. The line tension 
r is the excess grand-canonical free energy per unit length of the contact line, 

Lt = U £ -J2^n!. (3.7) 

i 

From this and from (3.6) for dU e , with Aj = N?/L, one obtains as the linear adsorption equation 

dr = — Aj dfii + c a da + cpd(5 + Cjd'y (3.8) 
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in place of (3.1). 



Another view of the physical origin of the additional terms c a da + • • • in (3.8), which makes it clear 
why the terms that must be added to — £ Aj dm are proportional to the infinitesimal changes in the contact 
angles that accompany the changes dm, is in the mean-field density-functional theory discussed in Section 
4, below. There it is also seen how in principle one may calculate dr + E Aj dm, i.e., c a da + cpdf3 + c^dj, 
from a microscopic theory. 

The invariance of r to the choice of location of the contact line noted in Section 2 now implies that 
the dependence of c a da + • • • on that choice must be just such as to cancel the contact-line dependence in 
£ Aj dm that was remarked in (3.3). Thus, if Ar is a shift in the location of the contact line and Ac a , Ac^, 
and Ac 7 are the resulting changes in the coefficients c a , etc., we must have 



Acq, da + Acg d(3 + Ac 7 dj 



■(e a p da af 3 + epy da^ + e« 7 da aj ) ■ Ar, 



(3.9) 



now with da etc. and da a p etc., as in (3.3) and (3.4), the infinitesimal changes in the equilibrium contact 
angles and interfacial tensions that accompany the infinitesimal changes dm- By (2.6), the right-hand side 
of (3.9) may equally well be written [o a p de a p + dep 1 + a ai de ai ) ■ Ar. 



Let us again fix fi4, . . . , m-+i, so that from (3.8), 



dr 

dm 



da 



dp dj 



- Al ^ + Ca diT l +c "diT l +c ^ 1 



(3.10) 



Thus, it is not — A 1( - 2 ,3) that is invariant to which of the two or four intersections of the hyperbolas A 2 = 
and A3 = is taken for the location of the contact line, but rather — A 1 ( 2 ,3) + c a da/ dm + • • • . This, with 
(3.9), implies that the combination 

Al < 2 < 3 ) + (^f ^ + e * + ^7 *~) • (r23 " ro) ' (3 - n) 

with ro any fixed location in the plane of Figs. 3 and 4, is independent of which of the intersections r 2 3 of 
A 2 = with A 3 = is taken. 



From (1.1), 
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rf(ro) 



rf (r ) a " 2 



d]*2 
dm 



rf (r ) " M3 



dm 

dm 



(3.12) 



and similarly for da ^ I dm an d da aj /dm', an d where ro, which is arbitrary because of the invariance of 
a a/ 3 to the choice of dividing surface, may be conveniently chosen to be the same as in (3.11) and to be 
the origin of the x, y coordinate system in the equations (2.7) for the hyperbolas Aj = 0. The quantities 
dmi 'dm an d dm/dm i n (3.12) are expressible in terms of the bulk-phase densities pf , , pj for i = 1, 2, 3 
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(3.14) 



9 



with 



A 



P2 -4 P3- P3 



7 P 7 
P2 ~ P2 P3 ~ P3 



(3.15) 



We may now illustrate the invariance of (3.11) to the choice of ^3; i.e., to the choice of which of the 
two or four intersections of the hyperbolas A2 = and A3 = one takes in (3.11) when da a p/dpi etc. 
and d^/dni and dps /dpi are as in (3. 12)-(3. 15). For this purpose we again adopt for the hyperbolas 
A2 = and A3 = those in Figs. 6 and 7 as determined by (2.10) and (2.11), the former for a case of two 
intersections and the latter for a case of four; and we again take the hyperbola in Figs. 8 and 9 as determined 
by (3.5) for our example of the hyperbola Ai = 0. Now from (2.12) and (3. 12)-(3. 15), and again with the 
numerical values quoted in the Appendix, we find for the expression (3.11) the common value —13.7 when 
r23 is identified as either of the two intersections PI or P2 in Figs. 6 and 8, again for every point ri on both 
branches of Ai = 0; and the common value —10.65 when r-23 is identified as any of the four intersections 
PI, ... , PA in Figs. 7 and 9, for every point ri on Ai = 0. Clearly, then, it is the full expression (3.11), but, 
as illustrated earlier, not A 1 ( 2 ,3) alone, that is invariant to which intersection ^3 of A2 = with A3 = is 
chosen as the location of the contact line in the definition of A 1 ( 2 ,3)- 



In the following Section 4 we view these matters in the context of a mean-field density-functional theory, 
as an example of a microscopic theory with which one may in principle calculate the coefficients c a etc. in 
(3.8). 



4 Mean-field density-functional theory 



In a general density-functional theory the line tension r would be obtained as 

/ & da - (a a pR a p + apyRpy + aajRt 
J A 



t = mm 



(4.1) 



where A is the area of the Neumann triangle in Fig. 4, da is an element of area, R a $ etc. are the lengths as 
labeled in that figure, and ^ is a functional of the spatially varying densities pi(r). The designation " min " 
means minimization over the pi(r), with r a vector in the plane of the figure. It is also to be understood that 
the limit in which the R a p etc. become infinite, and do so proportionally to each other, is taken in (4.1). 



In a general mean-field theory, and still with the convention that the entropy density s is one of the pi 
and the temperature T is the associated p, L , the functional ^ would be of the form 

* = F[pi(r),...,p c+1 (r); p u . . . , p c+1 ] + K, (4.2) 

where F is a local functional of the ft(r) for given pi,..., p c+ \ while K is a non-local functional of the 
Pi(r) but is independent of the pi. In this general mean-field theory F depends on the densities pi and the 
conjugate field variables p, L explicitly via 

F (pi, . . . ,p c+1 ; pi, . . .,p c+ i) = e (pi, . . . ,p c+ i) - ^2 PiPi+p{pi, ■ ■ ■ ,Mc+i), (4-3) 
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where e is the mean-field energy density as an analytic function of the densities pi, including the unstable 
parts, i.e., not reconstructed by the common-tangent-plane or convex-envelope construction; while p is the 
equilibrium pressure as a function of the pi. The terms — £ pipi +p in (4.3) correspond to the subtraction of 
that common tangent plane from e(pi, . . . , p c +i)- Thus, F is positive for all values of the pi except when, 
for the given pi, all have their bulk-phase values pf or p? or pj, when F, along with its first derivatives with 
respect to each of the p, L , are 0. 



Now let the thermodynamic state change infinitesimally via the infinitesimal changes dpi [only c — 1 
of which are independent because of the two Clapeyron equations (1.2) and (3.2)]. In then calculating the 
resulting dr from (4.1) we may ignore the contributions from the changes in the minimizing Pi(r) because 
r is stationary with respect to them; we need consider only the explicit dependence of the \£ on the pi 
contributed by the terms — £ pipi + p(pi, . . . , p c +i) in F. For the moment, we ignore the fact that the 
contact angles and so also the shape of the Neumann triangle change with the dpi, so that the region A of 
integration in (4.1) is for now taken to be unchanged. Taking account of 



dp = J2 Pi d ^ = Y^ p i d Pi = Y^ Pi dili 



(4.4) 



from the Gibbs-Duhem equation and the Clapeyron equations (1.2) and (3.2), and of the adsorption equation 
(1.1), we then have 



dT = YA\ [~^) + P?]da+ [ - Pi ( 
j I JA a JA L 



r)+p[ 



+ (TfR aP + TfRfr + T^R aj ) }dp t , 



da+ [-Pi(r) + pj] da 

J A~f 



(4.5) 



where A a , Ap, A 7 are the parts of the area A as so labeled in Fig. 4. But the linear adsorptions A« are 



A,; 



/ [Pi( 
JA a 



r)-pf]da + 



da+ j [pi(r) - p\\ da 

A-, 



Pi(r) ~ P 

-rfR ,-r/ /; -17 //,,.. 



(4.6) 



again in the implicit limit in which the R a p etc. all become infinite and do so proportionally to each other. 
Then (4.5) and (4.6) imply the (incorrect) linear adsorption equation (3.1). 



Here it is clear that what is so far missing from dr in (4.5) are the contributions from the changes in the 
areas of integration A a ,Ap, Aj that accompany the changes dpi in the thermodynamic state of the three- 
phase system, and which we have so far neglected. These changes in the areas result from changes in the 
contact angles, so the missing terms in (4.5) may be taken to be of the form c a da + cpdf3 + Cjd'y, 
as in (3.8). With a microscopic theory such as this mean-field density-functional theory, these additional 
terms may in principle be calculated because dr and EAj dpi are separately calculable. For example, with 
P4,..., p c+ \ again fixed, as in (3. 10)-(3. 15), c a da/dp\ + ■ ■ ■ may be evaluated from 

da dp dj dT dpi dp 3 

c a - Vcp- h c 7 — = h Ai + A 2 h A 3 - — (4.7) 

dpi dpi dpi dpi dpi dpi 

with dpi/dpi and dps/dpi obtained from (3.13)-(3.15). Such a calculation is now being undertaken [11]. 
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5 Summary and conclusions 



In any plane perpendicular to a three-phase contact line the locus of choices for the location of that line 
that makes the linear adsorption A, of any one of the components i vanish is a rectangular hyperbola given 
by (2.7) with (2.8). Two such hyperbolas Aj = and A^ = may intersect in two or in four points. 
The relative adsorption A^ k y which is the adsorption of component i when the contact line is so placed 
that the adsorptions Aj and Afc both vanish, i.e., at one of the two or four intersections of Aj = with 
Afc = 0, is given by (2.12). Unlike its analog da = — £ Tj dpi for an interface, the quantity — £ Aj dpi for 
a three-phase line is not invariant to the locations of the contact line and dividing surfaces, and so cannot 
be the differential dr of the necessarily invariant line tension r. This is illustrated in Section 3, where it is 
seen that A 1 ( 2 ,3) takes different values when calculated with the two or four different intersections of the 
hyperbolas A 2 = and A3 = 0. When the adsorption equation is supplemented with terms that cancel 
the non-invariant parts of S Aj dpi, in (3.8), the resulting expression containing A^^y in (3.11), takes the 
same values whichever of the intersections of A 2 = with A3 = is taken for the location of the contact 
line, thus verifying the correctness of the modified adsorption equation. The added terms vanish in the limit 
in which one of the phases wets the interface between the other two. 

The supplemental terms in the adsorption equation arise from the changes da, d(3, dj in the contact 
angles accompanying the changes dpi in the thermodynamic state. From a general mean-field density- 
functional formulation of the calculation of line tension and the associated linear adsorptions, in Section 
4, one sees explicitly why those added terms are connected with changes in the contact angles; i.e., with 
the dependence on thermodynamic state of the shape of the Neumann triangle in Fig. 4. Since both the 
line tension r and the linear adsorptions Aj are explicitly calculable with such a density-functional theory, 
it would in principle allow the supplemental terms, which are missing from the earlier formulation of the 
adsorption equation [8,9], to be calculated. 
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Appendix 

Here we provide values for the physical parameters a, 7, pf, T^(tq), Aj(ro), etc. on the right-hand sides 
of (2.8) that are consistent with the Ai/Fi, Bi/F i: Di/F u and Ei/Fi in (2.10), (2.11), and (3.5), which 
determine the illustrative hyperbolas in Figs. 6-9. The physical parameters are needed for the calculation of 
the values of A 1 ( 2 ,3) f rom (2.12) and the values of the expression (3.11), with (3.12), all as quoted in Section 
3. 

All parameters and the results of all calculations are given in dimensionless terms, so each pi is in units 
of the reciprocal of a microscopic volume, each a in units of a microscopic energy per unit area, r in units 
of a microscopic energy per unit length, each V in units of the reciprocal of a microscopic area, and each A 
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in units of the reciprocal of a microscopic length. 

The contact angles were chosen to be 

a = 7 = 0.7tt = 126°, (3 = 0.6tt = 108°. (A.l) 

Of the remaining parameters on the right-hand sides of (2.8), for each i six are independent: pf — 
Pi> Pi ~ PL r f( r o), rf 7 (r ), rf 7 (r ), and A;(r ). These are to be consistent, via (2.8), with the 
choices of the four ratios Ai/Fi, Bi/Fi, DijFi, and Ei/Fi specified in (2.10), (2.11), or (3.5), but two may 
still be specified at will. For all the calculations in the text we chose these to be 

Ai(ro)=r^(ro) = l (i = 1,2,3). (A.2) 

The values of the remaining parameters are now determined assuming the values already assigned to the 
Ai/Fi etc. in (2.10), (2.11), and (3.5), and are in Table 1. 

The coordinates x, y of the intersections PI and P2 in Fig. 6 and of the intersections PI, . . . , P4 in 
Fig. 7 are in Table 2. To verify convincingly that the determinant in (2.9) vanishes requires that PI, . . . , P4 
be known with great precision. With the precision to which they are given in Table 2, MAPLE evaluates the 
determinant as —2 x 10~ 18 . 
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Table 1: Numerical values of the parameters in (2.8) that are consistent with (2.10), (2.11), or (3.5), and 
(A.1) and (A.2). 







Fnr f 9 1 1 1 

1 Ul \Z»- 1 1 ) 


Fnr H 51 


a P 

Pi Pi 






1 06866 
i .uuouu 


(3 7 
Pi ~ Pi 






9 7S1 OO 


1 1 l r o) 








1 1 l r o) 






a xyfOA 


a (3 
P? - P? 


-2.19046 


-2.10292 




13 7 

P2 - Pi 


1.41539 


2.10292 




rf(ro) 


5.25671 


0.850651 




rf 7 (r ) 


1.5485 


0.850651 




Ps-Pi 


0.152786 


-0.057506 




(3 7 
P3 ~ P3 


0.152786 


0.363079 




rf(ro) 


0.850651 


0.247214 




rf(ro) 


0.850651 


-0.247214 
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Table 2: Coordinates x, y of the intersections PI, P2 marked in Fig. 6 and PI, P2, P3, PA marked in Fig. 7. 



For Fig. 6, 


PI 


a,, y 




-3 355 -1 490 




P2 


x,y 


= 


1.013, 4.937 


For Fig. 7, 


PI 


X 




-3.951 245 182 973 155 802 3 






y 




-4.075 823 658 595 715 972 8 




P2 


X 




1.157 707 701 169 996 880 






y 




1.529 799 699 747 754 165 2 




P3 


X 




0.570 879 622 026 697 738 57 






y 




-1.151 478 850 368 232 023 4 




PA 


X 




6.222 657 859 776 461 183 7 






y 




-6.302 497 190 783 806 169 
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Figure Captions 



Fig. 1. Phases a and (3 with the diffuse interface between them, and with a Gibbs dividing surface 
chosen at an arbitrary height zq. 

Fig. 2. Three phases a, (3, 7 and their interfaces a(3, (3j, and 07 meeting at a line of common contact. 
The phases lie in the angles between the planar interfaces. The dihedral angles between planes, also called 
a, (3, 7, are the contact angles. 

Fig. 3. Two alternative choices for the location of the three-phase contact line. The view is along that 
line, which is perpendicular to the plane of the figure and so appears as a point, while the interfaces appear 
as lines meeting at that point. 

Fig. 4. The Neumann triangle, with sides perpendicular to the interfaces, intersecting them at the 
distances R a p, R/3j, R a ^ from the contact line. The interfaces divide the triangle into the areas A a , Ap, Ay. 

Fig. 5. Unit vectors e a/ g etc. in the directions of the interfaces and unit vectors g a/ 3 etc. rotated 
counterclockwise from them by 90°. 

Fig. 6. Two rectangular hyperbolas making two intersections, marked as P1,P2. Their coordinates 
(xi,yi) and (x2, 2/2) are given in the Appendix. 

Fig. 7. Two rectangular hyperbolas making four intersections, marked as P1,P2,P3,P4. Their 
coordinates (x±, yi), . . . , (2:4, 1/4) are given in the Appendix. 

Fig. 8. Rectangular hyperbola Ai = 0, with points PI and P2 from Fig. 6. 

Fig. 9. The same rectangular hyperbola Ai = as in Fig. 8, with points PI, P2, P3, P4 from Fig. 7. 
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